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Abstract 

We consider an arbitrary atomic system (n-level atom or many such atoms) interacting 
with a strong resonant quantum field. The approximate evolution operator for a quantum 
field case can be produced from the atomic evolution operator in an external classical field by 
a ‘quantization prescription’, passing the operator arguments to Wigner D-functions. Many 
important phenomena arising from the quantum nature of the field can be described by such 
a way. 


1 Introduction 

The behaviour of atomic systems interacting with a quantized electromagnetic field in a cavity 
has been studied for a long time. Even the simplest model of a single two-level system interacting 
with a single field mode in a lossless cavity, the Jaynes-Cummings model (JCM) [1], reveals 
interesting properties like collapses and revivals of atomic inversion oscillations, trapping states, 
Schrodinger Cats, etc. [1-3]. They stem from the nonlinear nature of the JCM. 

For a linear system, the Hamiltonian is usually a linear function of generators of some rep- 
resentation of a finite-dimensional Lie algebra, and the Evolution operator (EO) belongs to the 
corresponding Lie group representation. This property is referred as a Dynamical Symmetry. The 
most familiar examples are the Harmonic oscillator and the spin rotation in an external magnetic 
field. Their dynamical groups are, correspondingly, SU( 1,1) and SU( 2). Dynamical symmetry 
results in an equidistant spectrum (or a spectrum consisted of several equidistant parts), coherent 
states (i.e., nonspreading wave packets) and many other attributes of ‘harmonic behaviour’ [4, 5]. 

Real systems are oftenly nonlinear. Interesting phenomena appear, if the nonlinearity is ‘weak’, 
i.e., the dynamics is ‘almost harmonic’ one. For instance, the JCM in the classic field limit 
is equivalent to the spin-1/2 rotation in the external field, that gives the simplest example of 
the dynamical symmetry (two-dimensional representation of SU( 2)). JCM collapses and revivals 
appear in the case of a strong quantum field, when the system is slightly nonlinear. 

JCM dynamics is very instructive and we may say more about it. The JCM possesses an 
exact solution and its EO can be found explicitly. It is 2 x 2-matrix with coefficients depending 
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on photon operators, see, e.g. [3]. On the other hand, the number of excitation N is a constant 
of motion in this model, and its EO can be written as a direct sum of infinitely many 2x2- 
matrices corresponding to different values of N. Every matrix is similar to the spin- 1/2 rotation 
around x-axis, but different blocks rotate with different frequencies, fi/v, called quantum Rabi 
frequencies. The spectrum consists from infinitely many pairs of levels separated by 20^, and 
anharmonicity results in a nonlinear dependence of quantum Rabi frequency on the excitation 
number, flw -— Q A^ (g is a coupling constant). Therefore, the JCM can be considered as an 
interesting example of a nonlinear system, howeverj>ossessing the exact dynamical symmetry. 

Treating the field classically, we put tt N -» 0 ~ gy/n, using the classical Rabi frequency 
rather than quantum ones. (Here n is the field intensity in units of photon number.) For 
a strong quantum field with Poisson photon distribution, different fiyv contribute to dynamics, 
n — < N < n + Vn. The difference in frequencies gy/^+An — gy/n^Kng ~ \fKn/n ~ g 

becomes important for times ~ g ~' . This is just the JCM cqllapse time. It is in the range of modern 
experimental possibilities, both for Rydberg atom micromasers and for optical microcavities (see 
the references in the review [1]). The frequency of revivals is proportional to g / \/fi. Thus, it is 
not surprising, that JCM collapses and revivals can be described in the frame of expansion over 
the inverse field intensity [3]. 

Natural generalizations of JCM involve more resonant levels and more atoms. Already a 
system of many two-level atoms (the Dicke model [6]) does not allow an exact solution. However, 
for very general class of atomic systems interacting with the quantum field under the Rotating 
Wave Approximation, the excitation number remains a constant of motion. If one of the field 
modes contains a lot of photons, we can neglect the other modes and develop a perturbation 
theory with the inverse excitation number as a small parameter. Precisely, the initial number of 
photons must be much larger then the maximum possible number of atomic excitations A. This 
program has been realized at the level of wave functions for the Dicke model in Refs. [7] and for 
more general systems in Ref. [8]. 

Here, we shall find an explicit form of the EO for an arbitrary atomic system interacting with a 
strong resonant quantum field. The atomic operators for such a system form a finite-dimensional 
representation of some compact Lie algebra (see, e.g., [9]). For the case of identical atoms, the 
algebra depends on the number of levels while the representation depends on the number of atoms 
and symmetry properties of initial atomic state under the atomic permutations. For instance, the 
case of A identical two-level atoms excited from a symmetric state leads to the (A+l)-dimensional 
representation of SU( 2). The exact EO can be written as a direct sum of finite dimensional-blocks 
with different excitation numbers. The dimension of blocks for the case N > A is determined 
by the atomic algebra representation, but the exact Hamiltonian in every block is a nonlinear 
function of the representation generators (see Eq. (9) below). However, we shall show that the 
dynamical symmetry can be restored in the zeroth- and the first-order approximations. It means 
that the evolution operator in every block can be well approximated by the operator from the 
corresponding Lie group representation. In the Dicke model case this approximate motion is a 
rotation of the collective atomic pseudospin (of the length A/2) around x-axis. Once again, the 
rotation frequencies in different blocks depend nonlinearly on the excitation number (see Eq. (10) 
below). Therefore, our zeroth-order approximation possesses a dynamical symmetry in the same 
sense as the exact JCM solution. The difference is that the motion in every block is described 
by the appropriate representation of atomic algebra instead of spinor SU( 2) representations for 
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JCM. Being restricted to JCM case our theory just reproduces its exact solution. 

The work is organized as follows. In the next section we shall describe the model. In Sec. 3 the 
asymptotic evolution operator for an arbitrary atomic system interacting with a strong quantum 
field is found. It has a matrix form in the atomic basis with coefficients depending on the photon 
operators. Let us remind, that the evolution operator for the atomic system in an external classical 
field ( semiclassical EO ) is a matrix of a finite rotation from the atomic group representation. The 
corresponding matrix elements are calculated by the group representation theory (see, e.g., [10]). 

Remarkably, the approximate ‘quantum’ EO can be obtained from the semiclassical one by a 
simple ‘quantization prescription’ (see the text before Eq. (23)), which introduces the operator 
arguments into the Wigner D-functions. Therefore, our results enable to write (without calcula- 
tions) approximate matrix elements of the quantum evolution operator as far as the ones of the 
semiclassical EO are known. 

We demonstrate the convenience of the form (23) for the EO in Sec. 4, where it is used to 
reproduce the wave functions found in Ref. [8] and to prove the approximate factorization of the 
system wave function for special initial conditions. Making use of Eq. (23) drastically simplifies the 
original proof of factorization [8]. The wave functions of Ref. [8] contain the information about 
collapses and revivals, trapping states and Schrodinger cats and provide the correct structure 
of the field quasiprobability distribution for the systems under study. Therefore, the proposed 
asymptotic form for the EO describes all these phenomena connected with the quantum nature 
of the field. 

2 Description of the model 

We shall work with the following Hamiltonian 

H = u){h + h) + V, V = g{aX + + a'X~). (1) 

Here, a,a\h = a f a are the photon annihilation, creation and number operators, describing a 
cavity quantized field mode with the frequency u>. g is the coupling constant with resonant atoms 
placed into the cavity, (we consider the exact resonance case for simplicity), h is the bare atomic 
system Hamiltonian determining the configuration of atomic levels, X± are atomic operators 
describing transitions between resonant levels and obeying commutation relations 

[h,X-] = -X-, [h,X + } = X + . (2) 

Hamiltonian (1) with conditions (2) describes quite general atomic system interacting with 
a resonant mode of quantized field under the Rotating Wave Approximation. The Dicke model 
corresponds to the particular case when an atomic system consists of A two-level atoms [6]. Then 
the operators h,X± belong to the (ATl)-dimensional representation of sr/(2) algebra and obey the 
additional commutation relation 

[X + ,X-}=2 h (3) 

The simplest case of the JCM corresponds to the two-dimensional representation of su(2). 

We define the basis of the atomic algebra representation as 

|fc,7 }at, 0 < k < A, A|fc,7) at = (k - C)\k,~f)aU ( 4 ) 
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where k is a number of excitations in the atomic system, A is its maximal value, C is a constant 
corresponding to the bottom energy level of the atomic system (fc = 0 ). 7 denotes all the other 

atomic indexes (atomic level populations), which are given by eigenvalues of the operators from 
the Cartan subalgebra. For example, in the case of su( 2 ) algebra we have the only Cartan operator 
h and there are no additional indexes. Then the maximum possible number of atomic excitations 
A is equal to the number of atoms A, and the bottom level of atomic energy is C = 4/2. In the 
presence of additional indexes (say, for three-level atoms) the bare atomic levels may be degenerate 
and the dimension of the representation is then larger than . 4 + 1 . 

It follows from Eqs. (2) that the excitation number operator commutes with the Hamiltonian 

N = a + a + h, [N,H} = 0. (5) 

It is useful to introduce the basis 

\N,k) = \N — k)f <g> \k) at , 0 <k<A, N\N,k) = (N-C)\N,k), (6) 

where | n)j is a Fock field states, | k) at is a bare atomic state (4). For a fixed value of N, the 
Hamiltonian ( 1 ) in the basis ( 6 ) is a finite dimensional matrix. It’s rows and columns are numerated 
by the indexes k and 7 . (In Eq. ( 6 ) and below, we omit the index 7 .) 

We shall explore the field phase operators defined as [11] 

a * = exp(— ifyy/h + l, a = \Zn+Texp 

exp(-i<j>)\n) j = |n+l)/ exp(i</>)|n)y = |n-l)/, n > 0 . ( 7 ) 

These phase operators are unitary and provide the correct physical results for large rc’s. The 
following commutation relations are valid: 

f(h) exp {-i4>) = exp (-z>)/(n + 1 ), exp(i^)/(n) = f(h + 1 ) exp(f^), (8) 

where /(n) is an arbitrary function of the photon number operator (determined by its Taylor 
expansion). Eqs. ( 8 ) follow from similar relations for the operators a, a T 


3 The Evolution Operator for the strong field case 

We consider here the case of initially strong field. Then, the total number of excitations is 
larger than the maximum possible atomic excitation number, i.e.“N ~ n h ” , and we can build 

a perturbation theory with N~ x as a small parameter. It is convenient to divide the derivation 
into subsequent steps. 

1. We eliminate a , a * and h from the Hamiltonian (1) using Eqs. (7) and the excitation number 
conservation, h = N — h, 

V = g h+l exp(?^)A r + + X_ exp(-i<^)V^V-^ + l.j (9) 

2. The square roots in the last equation can be written as series of powers of the operator 

(N + 1/2)— 1 = g 2 /(l 2 , n(N) = gjti + 1/2. (10) 
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We get 

v = n 

H 0 = exp X+ 4- exp (-*</>) X-, (12) 

ffi = -l/2{h,H 0 } + ,.... (13) 

Here {,}+ stands for anticommutator. 

Our choice of a small parameter (Eq. (10) rather than iV -1 ) is important, since it provides 
vanishing first order corrections to eigenfrequencies, as it has been shown in Ref. [7]. In fact, 
the first-order corrections are in such a way included into the zeroth-order ones, that improves 
the quality of the zeroth-order approximation. According to Eq. (11), the whole time scale is 
determined by the factor f 2 ( TV ) , which plays the role of the Rabi frequency for the problem under 
study. 

3. It is convenient to use the following transformation 

Q = exp ^i4>(h + C)] . (14) 

This operator is unitary on the states n > A > k. Acting on the basis vectors it gives 

Q\k) at ® \n)f = exp(i<i>k)\k) ai ® \n) } = \k) at ® |n - k) } = \N = n, k). (15) 

Since the operator h + C has an integer spectrum, the operator Q is a direct sum of different 
powers of the phase operator. 

Directly from the definitions (14), (7) and from the commutators (2), (8) we find 

f{n)Q~ l = Q-'f{h + h + C), Qf(n) = f(h + h + C)Q, (16) 

QX + Q~ l = exp(^)X+, QX-Q~ l = exp(-i<j>)X+. (17) 

4. Here we shall find the zeroth-order EO. (Wave functions in the second order for the Dicke 
model have been calculated in Ref. [7].) Applying the Q - transformation to the zeroth-order 
Hamiltonian and using Eqs. (16), (17) we have 

V ~ Cl(h+h) = Q Cl(h— C) H c iQ~ l , H c i = X+ + X-. (18) 

Therefore, (^-transformation removes the nondiagonal photon operators from the zeroth-order 
Hamiltonian. It transforms the operator H 0 into the purely atomic operator H ch which has a 
sense of the atomic system Hamiltonian in an external constant classical field. Moreover, the 
(^-transformation removes the field operators from the coefficients i/ 0 , Hi , ■ v in all the orders. Si- 
multaneously, the (^-transformation changes the expansion parameter, f l(N ) — * Cl(h-C), trans- 
forming it into the function of the photon number operator. Therefore, the (^-transformation 
separates field and atomic variables in the expansion (11) in such a way, that atomic operators ap- 
pear in the coefficients and the photon number operator is included into the expansion parameter. 
Thus, the calculation of high order corrections involve only the atomic operators. 

5. Now we can calculate the matrix elements of the EO with the Hamiltonian (18) 

exp(-zH>) w exp [-itQ(JV)tf 0 ] = Q exp [-ft(n - C)H cl ] Q~ x (19) 


Ho + + 


( 11 ) 
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between the atomic states (they are still operators in the field space). Taking into account the 
action of the (^-operator on the atomic states, we find 

at (k\ exp(—itV)\m) a t « e^ k at (k\ exp [—itfl(h - C)H cl j \m) at e~^ m 

— at(k I exp I m) at exp(z<^(A; - m)), ( 20 ) 

where we have denoted 

T k — C) = g \J fi T k C T 1/2. (21) 

Shifting e x< ^ k to the right in the last line of Eq. (20), we have used the commutator (7). 

What profit have we got with Eq. (20)? We have separated the field phase operators, writing 
them on the right. The operator FlkHci contains the only field operator ft which commutes with 
all the other ingredients and may be treated as a Q-number in the calculation of the exponent. 
Thus, we can reconstruct the quantum field EO if we know the atomic EO in the external classical 
field 

U d (u>) = exp (-ituH cl ). ( 22 ) 

The matrix U c i(lj) is the rotation operator from the atomic group representation. It’s matrix 
elements are known from the standard group representation theory. For the Dicke model case, the 
operator H c i = X+ + X _ — 2 S x is just a generator of rotations around x-axis, and the matrix 
elements of exp( — itu>2S x ) are the usual Wigner ZTfunctions. 

The quantum EO matrix elements can be produced from the matrix U c i(u;) by the following 
‘quantization prescription’: 

(i) to substitute the group parameter u> in the row of the matrix U c i{uo) by the 

operator u) — ► = yjh + k — (7 + 1 / 2 , (here k is the atomic excitation number in the row (^, 7 ) ); 

(ii) to multiply every matrix element at {&, 7 |t/|m, Y) a t from the right by the power of the 
phase operator exp(i<f)(k — m)). 

Writing explicitly in Eq. (20) the additional atomic indexes 7 , we have 

at(k)'y\U\m,'y f ) at = at(k,7\Uci{&k)\rn^') at exp (i<j>(k-m)). (23) 

This equation is our principal result. It corresponds to the wave functions found in Ref. [ 8 ]. Being 
restricted to the case of a single two-level atom, it gives the exact JCM EO (see, e.g., [3]). The 
quantization prescription formulated above can be easily generalized to include detuning [ 12 ], 

4 The wave function factorization 

The most remarkable feature of dynamics with the EO (23) is an approximate wave function 
factorization for special initial states. Let the field initially be in a coherent state, | in) j = |a), a = 
V^ne^, where n and 0 are the initial photon number and the phase of the field. (The number e ^ 
may not be confused with the phase operator exp(i<^)). The initial atomic state is taken as an 
eigenstate of the operator 

HM) = e‘*X+ + ei*X. = e'* A /r d e-‘« = A p |pW>„„ | M)«, = £>«.- 
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Here, H cl is determined by Eq. (18). We shall call A p and |p(^)) B< as semiclassical eigenvalues and 
eigenvectors. The vector \p) at is also an eigenvector of the semiclassical evolution operator (22) 

U c ,(u)\p)at = exp (~iuX p t) I p)af ( 24 ) 

It is known [13], that coherent states with large photon numbers are approximately eigenstates 
of the phase operators. More precisely, 

exp(ik<j>)\a) = exp(ik(j))\a) + 0(k/VH). (25) 

Now we act by the EO (23) to the initial state |in) = |a) ® \p) at ■ Multiplying the equation from 
the left by an arbitrary atomic vector (fc| ot and substituting the EO matrix element, Eq. (23) we 
have 

at(k\ U |a ) f <8>\p{4>)) at w J2 e iHk ~ C Kt(k\U c i(n k )\m) atat (m\p) at ®\a)j 

m 

= at {k\e'*( k - c) e-' a ^\p) a t®\a) f 

= at (k\exp ^-i<Apff\/n + fe + l/2^ l p(</>) )a<® |a)/- (26) 

Here, in the first line, the phase factor coming from the action of the field phase operator e'^ k ^ 
on the coherent state, Eq. (25) cancels the atomic phase factor a t(m.| p(<ft) ) e ^ 4 at{ m \p) at to 

give e'^ k ~ c ^ . In the second line we use the property (24), that is U c i{&k)\p) at = e >tXp ^ k \p)at- We 
stress, that the operator- valued group parameter H*, depends on the number of the row, so for 
rows with different k' s we must use different semiclassical EO’s with different values of the group 
parameter. Finally, in the last line of Eq. (26) we substitute e ^ ^ at{^\p)at at{k\ p(<fr ) ) at &nd 

a( {Jfc|e -,tAp ^* = a t(k \ exp( — itX p g\J n + ft + 1/2). Since a t{k\ is an arbitrary atomic state, we can 
write the wave function in the form 

U(t)\in) « exp -itX p gyfi + Ji + 1/2 |m). 


Introducing the notation fi = gxjn + 1/2 we can approximate the square root as 


y/n + U 1/2 « y/ft + 1/2 + , + o ■ 

2 Jn -F 1/2 

and we find the factorized wave function 


+ 1/2 



\m) = 

I *,(*)> = 
I A p (t)) = 


| * p (t)>® | A P (Q), 

exp (-itgX p y/h + l/2^ |a)/, 
exp — itg 2 X p h/(2Vt)^ \p) at , 


(27) 


This wave function 'has been found by a different way in Ref. [8] and used for the discussion 
of the trapping states, Collapses, revivals and Schrodinger Cat states for the systems under study. 

If different semiclassical eigenstates contribute to the initial state, the wave function is a 
superposition of corresponding factorized states (27). The evolution operator Eq. (23) describes 
the dynamics from any initial state, such that the initial photon number is much larger then the 
maximum number of atomic excitations. 
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